Extended Born— Infeld Dynamics and Cosmology 



M. Novello,^ M. Makler,!'^^^ l S. Werneck,'' and C. A. Romero^ 

^ Centra Brasileiro de Pesquisas Fisicas, 
Rua Xavier Sigaud, 150, CEP 22290-180, Rio de Janeiro, RJ, Brazil 
^ Observatorio Nacional - MCT, 
Rua Gal. Jose Cristino, 77, CEP 20921-400, Rio de Janeiro, RJ, Brazil 
^ Universidade Federal do Rio de Janeiro, Instituto de Fisica, 
CP. 68528, CEP 21941-972, Rio de Janeiro, RJ, Brazil 
'* Umversidade Federal do Rio de Janeiro, Observatorio do Valongo, 
Ladeira Pedro Antonio, 43, CEP 20080-090, Rio de Janeiro, RJ, Brazil 
^ Universidade Federal da Paraiba, Departamento de Fisica, 
CP. 5008, CEP 58051-970, Joao Pessoa, PB, Brazil 

We introduce an extension of the Born-Infeld action for a scalar field and show that it can act 
as unifying-dark-matter, providing an explanation for both structure formation and the accelerated 
expansion of the universe. We investigate the cosmological dynamics of this theory in a particular 
case, referred as the "Milne-Born-Infeld" (MBI) Lagrangian. We show that this model, whose 
equation of state has effectively a single free parameter, is consistent with recent type la supernovae 
data, providing a fit as good as for the ACDM model with the same number of degrees of freedom. 
Furthermore, this parameter is tightly constrained by current data, making the model easily testable 
with other observables. Contrary to previous candidates for unifying-dark-matter, the sound velocity 
of the MBI model is vanishing both close to the dark matter state as well as near the cosmological 
constant state. This could avoid the problems on the matter power spectrum that were present in 
previous adiabatic dark-matter/dark-energy unification models. We also present a short discussion 
on the causal propagation in nonlinear scalar field theories such as the one proposed here. 
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I. INTRODUCTION 

In the current standard cosmological model, it is often 
assumed that two unknown components dominate the 
dynamics of the universe: dark-matter and dark-energy. 
Dark-matter has been for several decades a fundamental 
ingredient in cosmological model building. Its need is rec- 
ognized, for instance, to explain galactic rotation curves, 
virial motions and X-ray gas emission in galaxy clusters, 
large scale flows, and gravitational lensing, all of which 
are consequences of its local clustering. This clustered 
component contributes to roughly a third of the critical 
density. 

On the turn to the century, compelling evidence 
appeared for a negative pressure component — dubbed 
dark-energy — capable of producing the accelerated cos- 
mic expansion implied by type la supernovae data (SNIa, 
see ref. 0) and providing around two thirds of the den- 
sity needed to explain the almost flatness indicated by 
cosmic microwave background radiation (CMBR) obser- 
vations 

Now, more evidence has been gathered to support the 
existence of such a negative pressure component. For ex- 
ample, new SNIa observations extending to higher red- 
shifts i3|, 1^ indicate that the universe was decelerating in 
the past, implying that the apparent acceleration cannot 
be due to simple systematic effects on the data. Inde- 
pendently, constraints on the age of the Universe Q also 
support the need of cosmic acceleration. The combina- 
tion of several cosmological observables, such as CMBR 



anisotropies jQ], the large-scale structure 0, and X-ray 
clusters |g, point to the same conclusion. 

These evidences were further strengthened by the re- 
ported detection of a late Integrated Sachs- Wolfe effect 
through the cross correlation of CMBR anisotropies and 
the galaxy distribution Q . In a nearly flat universe this 
effect is significative only if some matter component with 
non negligible relativistic pressure is present. 

Thus, assuming the validity of Einstein's Theory of 
Gravity, most of the pillars of modern cosmology seem to 
converge to a same conclusion: the universe is filled with 
a negative pressure component, besides dark-matter. It 
is nevertheless puzzling that the dynamics of our universe 
be dominated by two seemingly independent components 
whose nature is unknown. Although several candidates 
are available for the dark-matter (such as the neutralino 
and the axion fl^) and for the dark-ene rgy (such as the 
cosmological constant or a scalar field [ij), none has 
been detected in laboratory experiments. This situa- 
tion has motivated the search for unifying-dark-matter 
(UDM) models, where a single component is responsible 
for both structure formation and the current accelerated 
expansion of the universe. 

One example of UDM is given by the Chaplygin Gas 
and generalizations thereof Il3l Q, 0| , which were 
widely studied in the literature and tested against ob- 
servational data (see for instance and references 
therein). The Chaplygin gas can be obtained from the 
Born-Infeld action for a scalar field and has several phys- 
ical motivations. Nevertheless, the simplest implemen- 
tation of this model is in disagreement with large-scale 
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structure data 

In this work we introduce a new model for dark-energy 
and dark-matter unification based on an extension of the 
Born-Infeld Lagrangian. It is shown that this matter- 
energy component can drive the late time acceleration 
of the universe while behaving as dark-matter in early 
epochs and at present in high density regions. We focus 
on a subclass of such model which we refer as the "Milne- 
Born-Infeld" (MBI) Lagrangian. We study the dynamics 
of the MBI model and show that a de Sitter state is 
an attractor for the solutions of cosmological interest. 
Moreover, the "sound velocity" Cs — 5p/5p is close to 
zero both in the early universe and in the de Sitter epoch. 
This could solve the problem of instabilities present in 
previous unifying-dark-matter models [T^ . 

We show that cosmological observables in the MBI 
model, such as the luminosity distance, depend upon a 
single effective free parameter R (besides the "standard" 
cosmological parameters). We use recent data on 194 
SNIa to place a strong constraint on R. The MBI model 
provides as a good fit to the data as the ACDM model 
with the same number of free parameters (assuming a 
flat universe). 

The paper is organized as follows. We start by re- 
viewing the dynamics of purely kinetic noncanonical La- 
grangians for a scalar field and display its equivalence 
to a perfect fluid description (sec. We then discuss 
the effective geometry for nonlinear scalar field theories 
(sec. mill and review the Born-Infeld model in light of 
this approach (sec. IIV|I . In section we introduce an 
extension of the Born-Infeld Lagrangian and discuss its 
role as a unifying-dark-matter candidate. We turn to a 
particular case of this Lagrangian — the MBI model — 
and derive its corresponding equation of state in I VII The 
sound speed and equation of state parameter, as well as 
a comparison to other UDM models, are discussed in sec- 
tion The dynamics of homogeneous solutions of the 
MBI model is studied in section IVIIII In section Hxl it is 
explicitly shown that this model has a late time acceler- 
ating phase. The Hubble parameter as a function of the 
scale factor is derived in H and the results are applied 
to type la supernovae data in section IXII We sum up 
the results and present our concluding remarks in sec- 
tion IXllI Some speculations about the possible role of 
the proposed model in the early universe are discussed in 
the appendix. 



II. KINETIC LAGRANGIAN AND THE 
PERFECT FLUID REPRESENTATION 

In this work we shall consider scalar field models whose 
Lagrangian have a noncanonical form 



L = V{^)F{W)., 



(1) 



where W := W^ipW^ip. 

Noncanonical (i.e., non-quadratic) kinetic terms ap- 
pear rather frequently in low-energy effective La- 



grangians, generated by radiative corrections in the un- 
derlying true theory (see e.g. |0|). For example, in 
string and supergravity theories, nonlinear kinetic terms 
arise generically in the effective action describing mod- 
uli and massless degrees of freedom ,2^. Also, scalar 
fields are common many particle physics theories beyond 
the standard model. It is therefore not surprising that 
scalar fields with noncanonical kinetic terms could play 
a relevant role in an effective description of nature. If 
we expect the equations of motion in classical theories to 
be of second order, it is natural to consider only kinetic 
terms which are functions of the square of the gradient 
of the scalar field [23|. If it is further assumed that the 
Lagrangian is factorizable, such as in tachyon conden- 
sates 22], we are led to expression This type of La- 
grangian also appears as an effective theory for the low 
energy expansion of fluctuations around a ghost conden- 
sate (neglecting terms involving more than one derivative 
of (/?, see e.g. |23)- 

Lagrangians of the form were first introduced in 
the cosmological setting in the context of inflation 
l25l |. then as a dark-energy candidate |2^, and finally as 
UDM candidates, such as in the tachyon [23| and ghost 
condensate |23j | models for dark-matter and dark-energy 
unification.^ 

For simplicity, throughout this work we shall consider 
Lagrangians with constant potentials, i.e. of the type 



L = L{W). 



(2) 



Such models were the first ones investigated in the in- 
flationary setting [231 ■ To this family belongs the Born- 
Infeld action, which leads to the Chaplygin gas model. 
Some general properties of this type of Lag rangian in the 
context of UDM are discussed in refs. |2lLl29l|. 

Before closing this section let us display the equiva- 
lence of Lagrangian Q to a perfect fluid description. 
The standard definition of the energy-momentum tensor 
is given by 



T — 



(3) 



which, in the case of the scalar field with a purely kinetic 
Lagrangian ||2Jl, leads to 



= giLV + 2^^^ V^</? Xlyip , 



(4) 



where L^y dL/dW. As is well-known, any energy- 
momentum tensor can be written in terms of a fluid by 
the choice of a particular frame represented by an ob- 
server endowed with a four- velocity field v'^ , yielding the 
decomposition 



^ It is worth mentioning that other nonlinear theories, such as for 
spin 1 fields, may also play an important role in cosmology (see, 
e.g., ref. |28| and refs. therein). 
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where the ten independent quantities p,p, q" and tTq^ are 
obtained through the projections of T^i, onto v" and the 
space orthogonal to it, and iVap is traceless.^ 

In the frame commoving with the gradient of the field, 
defined by the normalized vector 

the energy-momentum tensor (0J is equivalent to that of 
a perfect fluid, since the heat flux and the anisotropic 
pressure tt^j, vanish identically. In this case, the non iden- 
tically zero quantities are the energy density and pres- 
sure, given by: 

p=-L + 2LwW, p = L. (5) 

Moreover, since both p and p are given as a function of 
W only, one of these relations can be inverted to provide 
p — p{p)- We see therefore that there can be no bulk 
viscosity. This shows that the dynamics given by any 
purely kinetic Lagrangian is equivalent to that of a perfect 
fluid, both in its energy-momentum tensor and equation 
of state (eos). 

III. CAUSALITY AND EFFECTIVE 
GEOMETRY 

Motivated by the recent interest on nonlinear scalar 
field theories in the cosmological context 0, 01 HE 
HE US HEl ) we find it useful to present a brief discussion 
on the issue of the causal propagation of the field. Our 
purpose is to study the behavior of discontinuities of the 
equations of motion around a fixed background solution. 
Instead of using the traditional perturbation method (the 
eikonal approximation), we shall use a more elegant for- 
malism synthesized in the work of Hadamard |30(| . In this 
method, the propagation of high-energy perturbations of 
the field is studied by following the evolution of the wave 
front, through which the field is continuous but its first 
derivative is not. To be specific, let a be the surface of 
discontinuity defined by the equation 

a{x^ ) = const . 

The discontinuity of a function J through the surface a 
will be represented by [J]a, and its definition is 

The discontinuities of the field and its first derivative are 
given by 



^ Explicitly, the scalars p and p (energy density and pressure) are 
defined by p := T°'^VctVp and p := —^h^^pT"'^, the heat flux is 
g° := h°'l^v~'Tpj, and the traceless symmetric anisotropic ten- 
sion is tt""'^ := h°'i^hl^''Tfj.^ + ph°'l^ , where = - v^v^. 



while for the second derivative we have 

where the vector is the normal to the surface of dis- 
continuity. Using these values in the equation of motion 
for the field (p, 

V^{LwV''p)^0, (6) 

we obtain 

kf, k^ g^^^. = , (7) 

where the effective metric is given by 

g^^f = Lw g^" + 2Lww V^^ V V (8) 

and g^'^ is the background metric. Only in the case of a 
linear theory L = W, the metric that controls the prop- 
agation of the discontinuities of the field coincides with 
the background metric. 

Therefore, the propagation of discontinuities of the 
scalar field, which we shall refer as scalaron, follows null 
curves in an effective metric that is not universal, but 
instead depends on the field configuration. We should 
emphasize that this property is quite generic for any non- 
linear field theory (see for instance refs. _31„ .32^, .3^ .3^ 
and references therein). 

In terms of the background geometry we can re-write 
the equation of propagation as 

k.k^g^'' = {k.V^^f ■ (9) 

This means that in the background geometry the 
scalaron behaves as time-like particles in cases in which 
LwLww < 0, and it behaves as tachyons in the cases in 
which LiyLww > 0- We shall came back to this issue in 
the appendix. 

IV. EXCEPTIONAL LAGRANGIAN: THE 
BORN-INFELD CASE 

The effective metric is an auxiliary instrument to sim- 
plify the description of the propagation of the disconti- 
nuities of nonlinear field theories. In nonlinear dynamics 
it does not correspond to the background geometry in 
which the field acts. From this property one could con- 
clude that such geometry should not appear in any aspect 
of the dynamics of the field. This is indeed the case in 
the most common theories, however, in certain nonlinear 
theories the dynamics is described in terms of quantities 
associated to the effective metric. For example, in Elec- 
trodynamics there is at least one exceptional case which 
is worth to point out: the nonlinear field theory proposed 
by Born and Infeld js^. In that theory, the Lagrangian 
is written in terms of the determinant of the effective 
metric L = (— det geffY^^ EE|- We shall show below 
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that a similar situation occurs for the nonhnear scalar 
field theory. In order to describe such special dynamics 
we must face the question: what is the property of the 
Lagrangian L{W) such that its corresponding dynamics 
is given only in terms of the effective metric? We shall 
follow the same lines as for the case of Electrodynamics 
[32]. 

We note that the (inverse) covariant metric ten- 
sor of the effective metric defined by the relation 
9t/j9eff = '^M' is given by 

eff 1 '2LwW „ V7 

ffpi ^-j—9^.^- J — 7j — N V^(^ . 

(10) 

A straightforward calculation shows that the evaluation 
of the determinant of the effective metric yields the re- 
sult: 

det 3^7^ = Ll, {Lw + 2W L^w) ■ (11) 

It then follows that the unique theory that can be writ- 
ten in terms of its associated effective metric is the one 
provided by the Born-Infeld like form: 

Lbi = -VbW + e, (12) 

where b and e are constants. Indeed, using H12|) into the 
expression of the determinant (|ll|l we obtain 

Lbi ^ const, {det g';//)^^\ (13) 

The energy density and pressure of the effective fluid 
description for the Born-Infeld dynamics (from eqs. IT^ 
and El 

p = e/\/bW + e, p= -VbW + e, (14) 

are such that the equation of state takes a very simple 
form'^ 

p = --. (15) 
P 

A fluid with the above eos is known as Chaplygin gas^ 

M 

The Born-Infeld action for a scalar field was first in- 
troduced by Heisenberg in 1939, in the context of me- 
son theory This model reappeared in the string 
theory setting, where the Born-Infeld Lagrangian is as- 
sociated to the parametrization invariant Nambu-Goto 
d-brane action in a (d + 1, 1) space-time The 



^ It is worth noting that, although Lagrangian 1121 has two free 
parameters, only one appears in the eos. 

* This name, which became popular in the d-brane setting, is off 
course inspired on the eos proposed by Chaplygin in 1904. In the 
original expression, p is only the internal energy, not the total 
relativistic one, and the weak energy condition is not violated, 
as opposed to what is needed in the cosmological case. 



Chaplygin gas has recently attracted much interest in 
the cosmolo gica l framework, both as a candidate for 
dark-energy |l2l liH l as well as for unifying-dark-matter 
|4^ . In fact, as can be seen from IjlSI) it has a 
negative pressure that could power the accelerated ex- 
pansion of the universe. Also, the pressure is negligi- 
ble for high density regions and "mechanical" pertur- 
bations are stable {dp/ dp = > 0). Therefore this 
fluid could also be responsible for structure formation. 
Because of these properties and its simplicity, equation 
(|15|l arises naturally as a toy model allowing dark-energy 
and dark-matter unification Ts^. The behavior of the 
Chaplygin gas, and simple generalizations thereof, as 
unifying-dark-matter (or simply quartessence^) was ex- 
tensively discussed in the literature (see, for instance, 
refs. [ll|ii,d5fcdallllil|4^|43fcd4i45|,,4a] and refer- 
ences therein). 

V. EXTENDED BORN-INFELD LAGRANGIAN 

As pointed out in the previous section, the Born-Infeld 
Lagrangian arises rather naturally in the context of a ge- 
ometric description for nonlinear scalar fields. Moreover, 
as mentioned above, besides its formal appeal, this the- 
ory has an interesting behavior in the cosmological con- 
text. In fact, eos 1)151) has provided the first model of 
dark-matter and dark-energy unification jl3l Il4j . Nev- 
ertheless, this model is ruled out by observations of 
weak gravitational lensing combined with the large-scale 
galaxy power spectrum |l3] and by the anisotropics on 
the CMBR Besides, it is also disfavored by a combi- 
nation of observational data on the background evolution 
of the universe 0, 0| . 

A generalization of the Born-Infeld Lagrangian was 
proposed in 113 and leads to the so-called Generalized 
Chaplygin gas |14| |. However, this model suffers from the 
same problems as the Chaplygin gas itself 0, ^3 ■ 

Nevertheless, the successful aspects of the Born-Infeld 
model as a candidate for dark-energy and dark-matter 
unification motivate the search for other models that still 
retain some of its simplicity at the same time avoiding 
its shortcomings. As shall be discussed in details in the 
forthcoming sections, such a model could be given by the 
Extended Born-Infeld (hereafter EBI) Lagrangian repre- 
sented by 

L = -^/aW^ + bW + e. (16) 

The original idea of Born and Infeld was to deal 
with a theory that describes the evolution of electromag- 
netic fields that are bounded from above. Here we ap- 



^ As in this scenario there is a single extra component besides 
"baryons", neutrinos, and radiation (as opposed to having 
both dark-matter and dark-energy), the unifying-dark-matter is 
dubbed as quartessence |41| . 
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ply such limitation idea to the case of a scalar field, as 
in Lagrangian (|12|l . However, instead of having only a 
maximum value we also require the existence of a mini- 
mum (such that a < and &^ > — 4ae). Therefore, W is 
defined between the two roots of the polynomial in (|16|l . 
From Lagrangian (|16|l . the energy density and pressure 
of the effective fluid are given by 



p=-L 
p = L. 



W{b + 2aW) 



(17) 
(18) 



The equation of motion Q derived from Lagrangian 
(Unj is 



(6 + 2aW) 
I 



0. 



(19) 



Therefore, a particular solution for the dynamics of the 
scalar field is given by a constant value of W , namely 



Wo 



b 

2a 



(20) 



In this case, the functional derivative of the Lagrangian 
Lw vanishes and consequently the tensor T^t^ becomes 
identical to the energy distribution of a cosmological con- 
stant. We shall call such a configuration the fundamental 
state. As we shall see later on, in the cosmological case 
there exists two de Sitter states associated with such con- 
stant Wq, one which is stable and the other unstable. 
Equation H19|l allows also the trivial solution if = const. 
which corresponds to another "vacuum" state. We shall 
see latter that, at least in the cosmological case, the phys- 
ically relevant vacuum is given by the fundamental state 
(HOI. 

Notice that, near the roots of the quadratic function, 
L ~ and the scalar field acts as pressureless matter. 
Besides, close to at least one of these roots the density is 
positive. Thus we have the possibility of a matter-energy 
component that interpolates between dust (close to a root 
of ll6|l and a cosmological constant like configuration (at 
the fundamental state). Moreover, the presence of L in 
the denominator in expression H17|) shows that the scalar 
field behaves as dust in high density regions. This situa- 
tion is typical from unifying-dark-matter models, where 
this component acts as dark-matter in clustered regions 
and as dark-energy in low density regions, providing and 
explanation for the accelerated expansion of the universe 
|13|. We shall investigate this model in the cosmological 
setting in the forthcoming sections. Our analysis will be 
focused on a special case of Lagrangian ((TB)l as described 
below. 



in which e = 0. In this situation one root corresponds to 
'A p — p — Q state and there is a positive root (we choose 
6 > 0) corresponding to a dust limit. As we shall see 
in section IVIIII this model allows a solution representing 
an expanding empty (p = 0) universe, which is known as 
the Milne solution. For this reason we refer to this model 
as the "Milne-Born-Infeld" (hereafter MBI) Lagrangian. 
In this case it is convenient to express (|16f) as 



(21) 



where A :— \/|a[ and := bA^"^. For this model, the 
whole admissible spectra of values for W is contained in 
the domain (0, S^). The hmdamental state is located at 
EV2. 

In the MBI model, the energy density and pressure 
(eqs. [TTlandlT^ are given by 



AW"^ 



P 



{Y? - W) 



= -A^w (s2-vr). 



(22) 
(23) 



Here the expression of the eos p (p) is rather complicated, 
as opposed to the Born-Infeld case (eq. [T3|l . It is more 
convenient to write p as a function of p: 

p = (ay? + V^^S4 - v)V4p- 

It is easy to see that, for high densities (p ^ ^S^), the 
equation of state goes asymptotically to the Chaplygin 
one 

A^Y^ 

p~ . 

P 

The density and pressure of the fundamental state are 
given by 



AY^ 

Pv ■= = -Pv ■ 



(24) 



In figure ^ we show p as a function of p in the interval 
Y^ > W > S^/2, from the dust to the fundamental state. 
As shall be shown in section IVIIII this is the relevant 
interval for cosmology. 

Before exploring the dynamics given by the MBI La- 
grangian (|^ . let us point out a fundamental difference 
with respect to the Born-Infeld case and the Extended- 
Born-Infeld with e ^ 0. In these two models, the La- 
grangian has a canonical form for small values of W 
(W < e/b and W < b/a). However, the MBI model 
does not have such a canonical limit. Nevertheless, we 
shall see that the MBI Lagrangian provides an interesting 
model for the dark sector of the universe. 



VI. THE MILNE-BORN-INFELD 
LAGRANGIAN 

In the reminder of this article we shall restrict our- 
selves to a subclass of Lagrangian Hlt)|) , namely the case 



VII. SOUND SPEED AND EQUATION OF 
STATE PARAMETER 

In the fluid interpretation, the important quantity that 
controls the level of stability of perturbations is the 
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FIG. 1; Equation of state of the MBI model. The dimension- 
less pressure p/(AE^), is plotted as a function of p/iyAY?') 
starting from the fundamental state pv = AT? /2. Also dis- 
played is the "vacuum" eos (dashed curve). 



"sound velocity" squared = {dp/dW) / {dp/dW) [H, 
which in the MBI model is given by 

In the interval Y? jl < T4^ < (from the dust to the de 
Sitter state) it is immediate to see that the quantity 
is always positive. Its maximum is given by f?^^^^-^ = 1/9 
and occurs at Wa — 3 the value at which the scalar 

field's self gravity starts being repulsive (see section llX)l . 
For W < < and the fluid is mechanically 

unstable. However, as we shall see in section IVlIII this 
is outside the physical domain and such values are never 
attained in the cosmological framework. 

In flgure|21we plot the sound velocity as a function of 
the density starting at the fundamental state H24|) . It is 
interesting to point out some features of the MBI model 
regarding the sound velocity. First, its maximum value 
is a third of the velocity of light, whereas in the Born- 
Infeld model it is equal to c. Most importantly, near 
the de Sitter state the sound velocity is zero. As shall 
be discussed later, our universe is currently approaching 
this de Sitter state, thus one might expect the present 
scalaron sound velocity to be small. This is in contrast 
to most quartessence models studied so far in the lit- 
erature, where the maximum sound velocity is attained 
precisely at the de Sitter state 0|. The relatively large 
value of the sound velocity at present produces strong 
oscillations and power suppression in the matter spec- 
trum, which render these models incompatible with cur- 



FIG. 2: The sound velocity squared as a function of the di- 
mensionless density p/ (AE^) in the MBI model. 

rent data^ [13, Ell. The MBI model could avoid this 
problem since the sound speed is significative only dur- 
ing some limited interval of densities, as shown in the 
figure. 

It is useful to see how the sound velocity evolves as 
a function of w := p/ p, which is sometimes know as the 
"equation of state parameter" . Both and w are the 
quantities that appear in the equations for the evolution 
of perturbations in the fluid representation (see e.g. ref. 
|48|). The Cg X ui relation allows also to make a compar- 
ison with other quartessence models. 

In the MBI case we have w = 1 — T,'^/W, such that 

^ _w{l + w) 



In the Born-Infeld model, we have ct = —w, while for the 
Generalized Chaplygin gas fl2. UM. flSl l4ll | the relation is 
cj. — —aw. Recently, another generalization of the Born- 
Infeld action was proposed, the n*'* order Dirac-Born- 
Injeld theory, given by L oc —(1 — W)^/" in which 
case Cg — —2w/ {n+ (n — 2)(1 + 2w)). Other two models 
of quartessence were investigated in [l^ l45l | , where the 
equations of state are given by p = —M^ exp (— ap/M*) 
andp = -Af4/ [log (p/M*)]". Here we assume adiabatic 
perturbations of these models such that = dp/dp. 



^ A part icular kind of entropy perturbations solves this problem 
Il8ll44l but, as discussed in section|n] models with purely kinetic 
Lagrangians for scalar fields are adiabatic. Even other scalar 
field models which have the same background evolution as the 
Chaplygin gas seem to be highly disfavored |4^ . 
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(a) -1 -0.8 -0.6 -0.4 -0.2 (b) -1 -0.8 -0.6 -0.4 -0.2 

w w 




FIG. 3: The sound velocity as a function of the equation of state parameter w for the MBI case and other quartessence 
models, (a) Results for the MBI model (thick line) and the Generalized Chaplygin gas, with a — 0.1,0.25,0.5, and 0.75 (from 
bottom to top, short dashed curves), and a — 1 (equivalent to the Born-Infeld model, long dashed line), (b) Exponential 
quartessence with a — 0.1, 0.5, 1, 2, and e (from bottom to top), (c) Logarithmic quartessence with a — 10"'^, 10"^, 10"^, and 
(from bottom to top), (d) Results for the n*'' order Dirac-Born-Infeld theory with n = 2, 5, 20, 100, and 500 (from top to 
bottom) . 



In figure|31we plot as a function of w for the models 
referred above. In the upper left panel we show the re- 
sults for the MBI (thick line), Born-Infeld (long dashed 
line), and the Generalized Chaplygin gas models (short 
dashed curves). In the upper right, the results for the 
exponential quartessence are displayed. In the lower left 
panel we plot curves for the logarithmic quartessence 
model. Finally, in the lower right panel of this figure 
we show curves for the n*'' order Dirac-Born-Infeld the- 
ory. From these figures, it is clear that all the models 
start with = at w = 0, where quartessence acts as 
CDM. Nevertheless, in the MBI case, the sound veloc- 
ity vanishes again when w = —1, whereas in the other 
models displayed in the figure the sound velocity attains 
a maximum at that point. This leads to oscillations and 
suppression in the power spectrum of density perturba- 



tions.^ Thus we expect that this issue of the power spec- 
trum will be alleviated in the MBI model. 



VIII. DYNAMICS OF SPATIALLY 
HOMOGENEOUS SOLUTIONS 



Now, let us focus on the dynamics of the MBI model 
for spatially homogeneous and isotropic configurations. 



^ This was shown explicitly for the Generalized Cha plygin Gas llTl 
and the Exponential and Logarithmic models llSl . The behavior 
of perturbations was not yet studied in the n*'' order Dirac- 
Born-Infeld theory l4il . 
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In this case the Raychaudhuri equation is given by* 

e + \e^ + \{p + ip)^Q. (26) 

the energy conservation equation is 

p+(p + p)0 = O, (27) 

and the expansion parameter can be written as = 3a/a, 
where a {t) is the scale factor. 

Combining equations and lP7|l we obtain a dy- 
namical system: 

f=F(p,0), ^ = G(p,0), (28) 

with 

^W{T?~wy 3 2 (S2 - W) ' 

(29) 

where, from eq. (|22|l . is given by 

2-i/3(3^SV)'^fl + Vl + 27fey^' - (2/3)^/^^ 
W= ^ ^ 

(9^4s2p2)l/3 

(30) 

From an inspection of H2t)l I27|) it is easy to find the four 
stationary points P(p, &) of this system. Two correspond 
to the fundamental state W = S^/2, 

M = {p^,^/3^), Ar = (pv,-v/3^), (31) 

and two are on the line — 0, 

O = (0,0), Q = (pc,0), (32) 

where pv is given by and pc — (3\/3/2)pv Consider- 
ation of the behavior of the system in the neighborhood 
of these points shows that M and N are two-tangent 
nodes (one stable and the other unstable) while those on 
the 9 = line are saddle points. 

In figure 0] we display the phase diagram in the {p, 9) 
plane for equations (|28I3U|) . as well as some trajectories. 
From expression H22I) we see that only positive energy- 
density solutions are allowed. From the figure, it is clear 
that any trajectory that starts with a density lower than 
the fundamental state will remain with p < Pv Con- 
versely, trajectories starting at p > pv are bounded to 
have p > Pv It is easy to see from the Friedmann equa- 
tion that the solutions for the fiat case {k = 0) give 
p = 9^/3, which is the dashed trajectory plotted in this 



For simplicity, in the following sections we shall consider a uni- 
verse composed only of scalaron matter (baryons will be included 
in sections lYl and [Xil . We choose units such that 8nG = 1. 



P 




FIG. 4: Vector field and trajectories of the MBI model in 
the (p, 9) plane. The dashed curve is the trajectory for a fiat 
universe. The dots correspond to the four stationary points 
M, A'', Q, and O discussed in the text. 



figure. Notice the existence of a trajectory with vanish- 
ing density and nonzero expansion (the bottom horizon- 
tal line). This is the well known Milne solution, and is 
allowed only if we have e = in Lagrangian , leading 
to our nomenclature for this model. 

From now on, we shall be concerned with the right 
hand side of figure 01 i.e., with expanding solutions. 
Some of these may correspond to the background (unper- 
turbed) cosmological solution, in particular those start- 
ing from high values of p. As a nearly flat spatial ge- 
ometry is favored by current astrophysical observations, 
one should only consider trajectories near the fc = one. 
From figure 21 it is clear that these trajectories end up in 
the fundamental state M. In other words, the de Sitter 
solution is an attractor for the trajectories of cosmologi- 
cal interest. 

It is important to notice that, even in the inhomo- 
geneous case, there is a "vacuum barrier" that cannot 
be transposed in any expanding region. This can be seen 
from the energy conservation equation (|27|l . which is valid 
for any 9. As the vacuum state p = —p is approached, 
the density freezes out, and this is a stable point (the 
contrary occurs in a contracting region). In the MBI 
model, only states that initially violate the weak energy 
condition are able to approach the point O, which rep- 
resents a "static vacuum" ((p = const.). Therefore, as 
pointed out before, only the domain between the dust 
regime W ^ and the fundamental state Wq = E^/2 
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— which represents a "dynamical vacuum" — is relevant 
for the cosmological model. We shall see below that, 
within this interval, the universe has an accelerated ex- 
pansion phase. 



IX. THE ACCELERATION OF THE UNIVERSE 

The acceleration of the universe is controlled not only 
by the energy density but also by the pressure, according 
to the Raychaudhuri equation (|26|l . Writing this equation 
in terms of the scale factor a and using equations (|22|l and 
(ESJ yields: 



AW {AW - 3S2) 
^JW{Y? - W) 



(33) 



Starting with W close to S^, that is, a matter dominated 
universe, it is clear that the acceleration changes sign at 



Wa^-T.' 



(34) 



If the scalaron is the dominant component of the cos- 
mological energy density, the universe will begin to ac- 
celerate when W is close to Wa- Therefore, along the 
transition from a dust universe to a de Sitter one, the 
expansion starts accelerating, as expected. Close to this 
transition point the sound velocity reaches its maximum 



value Cg — 1/9, as discussed in section IViT 



X. THE HUBBLE PARAMETER 

As discussed in the preceding sections, the MBI La- 
grangian appears as a interesting model for the matter 
content of the universe, providing an explanation for the 
accelerated expansion as well as for structure formation. 
Clearly, in the cosmological setting, among all solutions 
displayed in the phase plane of figure ^ we must con- 
centrate in the domain characterized by the inequality 
T?/2 <W <T.'^. 

To compare the predictions of the model to observa- 
tional data on the background dynamics, such as the 
redshift-distance relation, we need to know the evolu- 
tion of the Hubble parameter. For this sake, the energy 
density has to be expressed in terms of the scale factor. 
This is accomplished by solving the energy conservation 
equation (|27|1 . which is conveniently expressed in terms 
of W . Using (|22|l and after a few manipulations, we get 



1 



w 



2W 



Integrating the above equation, in the desired range of 
W , we obtain 



SVSCa 



2W ~T? 



(35) 



where C is a (conveniently defined) integration constant. 
This relation is easily inverted to get 



W = 



(36) 



were the sign of the root was chosen such that W 

when a — s- 0, to have a dust behavior in the young 
universe. Inserting H36(l in H22|l and using l|35|l . we finally 
get 



AY.' 



2 {\^Ca-^{^fYTC- 



1)] 



3/2 



1 



(37) 

,-3 



It is easy to see that, in the limit a^O, p—^ALZ?\J%Ca 
and the scalaron density scales as CDM. On the other 
hand, in the limit a— >cxd, the fundamental state (|24|l is 
reached. Thus, as expected, this model behaves as CDM 
at early times and goes to a cosmological constant-like 
state in the future, being a candidate for quartessence. 
In fact, this asymptotic behavior occurs in any equation 
of state that intersects the p = —p line [T^ . 

In a realistic cosmological situation one has to consider 
also the baryons and the radiation components. Defin- 
ing fij = Pi/pcrit = ^T^Gpil (3i/^), with H = a/a, the 
Friedmann equation is given by 



SttG 



nla- 



where = —k/HQ, the subscript i denotes each compo- 
nent, and the knot means that the quantity is computed 
at the present value of the scale factor a — ag 1. For 
baryons Pb/Ph = o,~^, while for radiation Pr/Pr — a~^- 
The energy density in the scalaron field (|37|l shall be de- 
noted as ps- From expression (|37|1 it is clear that Ps/ p^s 
a function only of C and a. Therefore, as in the BI case, 
although the Lagrangian H21|) has two free parameters, 
the observables depend upon a single free parameter. As 
the parameter C lacks a direct physical interpretation, it 
is convenient to define the quantity 



R 



A^ 



(38) 



which, as can be seen from (|37|) . depends only on C. The 
parameter i? is a measure of how close the present energy 
density in the scalaron field is to its final "vacuum" value. 
In other words it is a measure of how close our universe 
is to its final de Sitter state. From the definition above, 
it is clear that < i? < 1. The explicit form of C (i?) is 



1_(2/-1)^ 



1-/ 



with 



(3V3i? + Vl + 27fl2)'/'^l V3 
i?(3\/3i? + Vl + 27i?2)'/^ 6 



/ = 



(39) 



(40) 
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Defining Ps/ p1 ■= g (a, R), were g is obtained from 13 7() 
together with H39|) and (|40|l . the Hubble function will be 
given by 



0„-3 



(41) 



where 17° = 1 - (172 



f7g 



n^) and i/o is the Hubble 
constant. In the next section we shall use recent data 
on the distance-redshift relation determined from type 
la supernova observations to set constraints on R. 



XI. CONSTRAINTS FROM TYPE lA 
SUPERNOVAE 



Type la supernovae (SNIa) observations have now be- 
come a standard tool for cosmology. Empirical studies 
show that these supernovae can be good standard candles 
after hght curve calibration (see e.g. refs [13, E3)- The 
first compelling evidence from a single experiment for the 
accelerated expansion of the universe was indeed based 
on SNIa data . Larger data sets including higher red- 
shift supernova are currently available 3^ J., _52J . These 
new data makes a stronger case for the current accel- 
erated expansion, since, as expected, a deceleration is 
detected for higher redshift supernovae, weakening the 
concerns that the reported acceleration could be due to 
systematic eff'ects. 

Here we shall use some of these recent data to check 
whether the MBI model can provide a good fit to SNIa 
observations and to set constraints on the parameter R 
(eq. I38|l . For this sake, a sample of 253 SNIa was com- 
piled from references Q and [52] which provide tables 
with the redshift z, logi?L, and its variance crj^ogDi ^'^^ 
each supernovae. Here -Dl is the luminosity distance 
times the Hubble constant, which, in the fiat case, is 
given by 



dz' 

DL=dLHo = {l + z)H„ / —— 

Jo H{z') 



(42) 



The Hubble parameter H (z) is obtained from equation 
(eq. with z = - 1. 

Since observations of anisotropics in the CMBR indi- 
cate that the Universe is nearly flat (Sl^, ^ 1), the anal- 
ysis that follows will be restricted to the zero curvature 
case. Also, at the redshifts probed by SNIa data, the 
radiation contribution to the Hubble parameter (eq. I1T)| 
is negligible. Therefore, we set fi^ -I- fi^ = 1. We also fix 
the baryon density parameter at fif, = 0.04, in agreement 
with the observed abundances of light elements [sj Is^ 
together the Hubble Space Telescope key project result 
^5]. After these choices, Dl becomes solely a function 
of the MBI parameter R. 

Following and [H^ we discard local supernovae with 
z < 0.01, because the peculiar motion contribution to z is 
too high, and those with high host extinction. Ay > 0.5, 
which could cause a strong bias in the determination 



of D^. After these cuts, we end up with a sample of 
194 SNIa extending up to z = 1.75. Despite eliminating 
low redshift SNIa, peculiar motions affect the measure- 
ment of z at all redshifts, causing a scatter around the 
value given by the bulk cosmic expansion. This is taken 
into account by including an uncertainty in z which is 
propagated into the luminosity distance and is summed 
in quadrature with the observational uncertainty in D^. 
Therefore, the Chi-squared is defined by: 

^ g [log(jg^- (zQ ) - log {D^ {^^))] ' 



i=l 2 

^log(Z5i(z,)) 



* log D 



Th 



where the theoretical prediction is computed from 1)42(1 
together with (|41ll and the observational values are given 
in the tables of 3] and 52]. Here we assume a velocity 
dispersion of — bOOKm/s (such that — o^jc ~ 
1.7 X 10~5). 

With the above described assumptions and data set, 
the best fit value is i? = 0.713 which corresponds to = 
198.0. Therefore the Chi-squared per effective degrees 
of freedom (193 in this case) is given by Xdo/ = 1.026. 
The MBI model provides a fit as good as for both the 
ACDM (Xrfo/ = 1-029) and Born-Infeld (x^/ = 1-020) 
models, which also have a single free parameter. The 
reduced Chi-squared is also very close to other mod- 
els of quartessence with two free parameters, such as 
the Generalized Chaplygin Gas (x^^y = 1.025), Expo- 
nential Quartessence (Xdo/ = 1.021), and Logarithmic 
Quartessence (xL/ = 1-023) However, as pointed 

out before, these models are incompatible with large- 
scale structure data for adiabatic perturbations 0, 0| . 

Defining the likelihood by C.{R) oc we have 

checked that the probability distribution of R is very 
close to a Gaussian, with less than 1% deviation. At 
Icr (68.3%) confidence level we have R = 0.713 ± 0.029. 
Thus, the current SNIa data is able to limit the MBI free 
parameter to a rather narrow interval. This makes the 
model readily testable by using different observables. If 
the values of R are incompatible to the ones obtained 
from SNIa, the model must be discarded. 



XII. CONCLUDING REMARKS 

Inspired in the Born-Infcld theory for a scalar field, 
we have introduced a new model for dark-matter and 
dark-energy unification, the EBI Lagrangian H16|l . We 
have developed a systematic study of a subclass of this 
Lagrangian, namely the MBI model, for which the tra- 
jectories in the phase space of homogeneous cosmologies 
were computed. It was shown that the trajectories of cos- 
mological relevance end up asymptotically in a de Sitter 
state. 

A new feature of this model is the fact that the equa- 
tion of state changes concavity. Therefore, the sound 
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velocity is maximal at an earlier epoch, contrary to the 
previously considered convex quartessence models, where 
Cs is greater at present 0. This could solve the known 
problems associated with the power spectrum of adia- 
batic quartessence |0. Therefore, this model has ad- 
vantages over other generalizations of the Born-Infeld 
Lagrangian [iRIi^. 

We have obtained an explicit form for the Hubble func- 
tion, which depends on a single free parameter (once the 
fi^ are fixed), expressed in terms of the ratio of the funda- 
mental state and current scalaron densities {R = p^/ pi). 
We applied this result to obtain the luminosity distance- 
redshift relation and used recent type la supernovae data 
to place constraints on R assuming a flat universe. The 
model provides a good fit to the data, giving a similar 
Xdo/ ACDM model — with the same number of 

degrees of freedom — and other models of quartessence. 

We find that R = 0.713 ± 0.029, at 68.3% confidence. 
Therefore, contrary to other quartessence models , w hose 
allowed parameter values are highly degenerate [la, EE 
0, 1^ , this model is tightly constrained by current data. 
Once R is fixed, the model provides clear predictions for 
other cosmological observables, such as the matter and 
CMBR power spectrum. 

The results reported in this paper motivate further 
studies of the EBI model, both to perform a detailed 
study of its dynamics in the general case, as well as to 
obtain observational constraints on the model. In partic- 
ular, in what regards the MBI case, a more complete ob- 
servational analysis should be performed, both to check 
if the model withstands to other observables and to im- 
prove the statistical analysis. For example, we could 
have marginalized over the Hubble constant and baryon 
fraction, instead of fixing these parameters to their best- 
fitting values. We could also study the effects of spacial 
curvature on the determination of the parameters. It is 
fundamental to investigate the formation of cosmic struc- 
tures in this scenario, both in the linear and nonlinear 
regimes. Nonetheless, the results of this work indicate 
that the model here introduced may be a viable candi- 
date for quartessence. 

As a final remark, let us recall that in the quartessence 
model there is an effective transformation of dark-matter 
into dark-energy and that, in the homogeneous case, 
around 96% of the cosmic matter ends up in the de Sitter 
fundamental state. This brings us to a comment by S. 
Carrol, stating that "[the de Sitter solution] represents 
the only natural stable solution for cosmology, and one 
of the outstanding pro blems is why we don't find our- 
selves living there" (53 ■ In the MBI model, and in the 
quartessence scenario in general, the answer to this ques- 
tion is: we are indeed heading there! 
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APPENDIX: SPECULATIONS ON THE 
POSSIBLE ROLE OF SCALARONS IN THE 
EARLY UNIVERSE 

In this section, we conjecture about the possible role 
of scalarons in solving two fundamental problems of cos- 
mology: the homogenization of the universe and the issue 
of initial conditions. First, let us point out that scalarons 
propagate as tachyons. Indeed, from section llTTl the prop- 
agation of the scalarons is provided by k^kj^g^'^^ — 0, 
where the effective metric is given by equation Thus, 
in the domain T."^ > W > /2 we find, from (O, that 
the scalarons propagate outside the local light cone 

For the present energies in the cosmos, the interaction 
of the scalaron with normal matter must be very weak, 
otherwise it would not be dark or it would have already 
been detected in the laboratory. However, the coupling 
of the scalaron to ordinary matter could be important in 
the early universe. One may wonder if the super luminal 
propagation of scalarons could have an important role in 
the homogenization of the universe. 

Another interesting point relates to the initial condi- 
tions. One of the possible natural initial configurations 
is that the scalar field be produced in the state of maxi- 
mal "kinetic term", i.e., with maximal space-time fluctu- 
ations. If on one hand it may be hard to understand how 
could the universe be "prepared" in a state of infinite 
density, on the other hand it is perhaps natural to accept 
that the scalaron was "prepared" in such a state of max- 
imal kinetic term V^t/sV^tp = S^, which in turn leads 
to a divergent p. If radiation is coupled to the scalaron 
in the early universe, the temperature could be very high 
due to the thermal equilibrium with the scalaron. A cou- 
pling to the electromagnetic field could be introduced in 
a simple way considering a gauge interaction with a com- 
plex scalar field. Notice that the results of the previous 
sections do not change if we consider a complex scalar 
field instead of a real one. For a discussion of couplings 
to other fields see e.g. refs. [23ll58l| . 

Therefore, the MBI model might provide some hints 
about the early universe, being responsible for the ex- 
tremely high densities as well as for the cosmic homog- 
enization. These speculations will be subject of further 
study. 
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